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ABSTRACT 

The observational data on the anisotropy of the cosmic microwave background con¬ 
straints the scalar spectral tilt n s and the tensor to scalar ratio r which depend on the 
first and second derivatives of the inflaton potential. The information can be used to 
reconstruct the inflaton potential in the polynomial form up to some orders. However, 
for some classes of potentials, n s and r behave as n s (N ) and r(N) universally in terms 
of the number of e-folds N. The universal behaviour of n s (N) can be used to recon¬ 
struct a class of inflaton potentials. By parametrizing one of the parameters n s (N ), 
e(N) and 4>(N), and fitting the parameters in the models to the observational data, 
we obtain the constraints on the parameters and reconstruct the classes of the infla¬ 
tionary models which include the chaotic inflation, T-model, hilltop inflation, s-dual 
inflation, natural inflation and R 2 inflation. 
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1 INTRODUCTION 

The quantum fluctuation during inflation seeds the large- 
scale structure and imprints the information of early Uni¬ 
verse in the anisotropies of the cosmic microwave back¬ 
ground (CMBR). For the slow-roll inflation with a sin¬ 
gle scalar field, the scalar (Mukhanov & Chibisov 1981) 
and tensor (Starobinsky 1979) power spectra can be 
parametrized by the slow-roll parameters (Stewart & Lyth 
1993). In particular, the scalar spectral tilt n s , the ampli¬ 
tude of the scalar spectrum A s and the tensor to scalar 
ratio r can be constrained by the observational data on 
CMBR. For an inflationary model, we can calculate the 
slow-roll parameters and the observables n s and r, and 
compare the results with the observations. However, there 
are lots of inflationary models with different potentials 
(Martin et al. 2014a) and it is not an easy task to com¬ 
pare all the models and constrain the model parame¬ 
ters in the models, although we may use Bayesian ev¬ 
idence to select best models (Martin et al. 2014b). Re¬ 
cently, the temperature and polarization measurements on 
the CMBR by the Planck survey gave the results n s = 
0.9645±0.0049 (68% CL), ln(10 lo A s ) = 3.094T0.034 (68% 
CL) and ro .002 < 0.11 (95% CL) (Planck Collaboration I 
2015; Planck Collaboration XX 2015). If we take the num¬ 
ber of e-folds before the end of inflation at the hori¬ 
zon exit N = 60, then the measured scalar spectral 
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tilt can be approximated as n s — 1 ~ —2/A, and r 
can be described as r ~ N~ p with p > 1. The small 
value of r alleviates the problem imposed by the Lyth 
bound (Choudhury & Mazumdar 2014b; Gao et al. 2014; 
Choudhury 2015; Gao et al. 2015). Since the running of the 
scalar spectral index is a second-order effect, it is expected 
to be in the order of 1CU 3 for slow-roll inflationary models 
(Gao & Gong 2014). The running of the scalar spectral in¬ 
dex is constrained to be dn a /dlnk = —0.0057±0.0071 (68% 
CL) (Planck Collaboration XX 2015). 

For chaotic inflation with the power-law potential 
V((j>) = Vo<t> p (Linde 1983), we get n a — 1 = — (p + 2)/(2A) 
with N being the number of e-folds before the end of in¬ 
flation at the horizon exit. For the hilltop models with 
the potential V{(j>) = Vo(l — ^/> p ) (Boubekeur & Lyth 
2005), we get n a — 1 = —2 (p — l)/[(p — 2)A] and r ~ 
0. For the Starobinsky model (Starobinsky 1980), we get 
n a — 1 = —2/A and r = 12/A 2 . For a class of in¬ 
flationary models with non-minimal coupling to gravity 
(Kallosh et al. 2014), it was found that n s — 1 = —2/A 
and r = 12/A 2 . For natural inflation (Freese et al. 1990), 
we get n a - 1 = -[1 ± exp(-A// 2 )]/[/ 2 - f 2 exp(-A// 2 )] 
and r = 8exp(-A// 2 )/[/ 2 - f 2 exp(-A// 2 )]. Therefore, 
there are some universal behaviours for n a and r which 
are consistent with the observations for a large class of in¬ 
flationary models. The suggested simple relation between 
n a and A tells us that it will be more easier to com¬ 
pare the models with the observational data if the slow- 
roll parameters and the observable n a can be expressed 
in terms of A. By parametrizing the slow-roll parameter 
e as t = /3/(A ± 1)“, Mukhanov derived the correspond- 
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ing infiaton potential (Mukhanov 2013). For the simple in¬ 
verse power-law form N~ p , Roset divided some inflation¬ 
ary models into two universal classes with the same 1/N 
behaviour for n 3 and different power-law behaviour for r 
(Roest 2014). More complicated forms of e(N) were also pro¬ 
posed by Garcia-Bellido and Roset (Garcia-Bellido & Roest 
2014). The parameters a and /3 in the parametrization were 
fitted to the observational data in (Barranco et al. 2014). 
The consequences of the two fixed points on n B and r in the 
parametrization were discussed in (Boubekeur et al. 2015). 
In addition to the parametrization of e, the reconstruc¬ 
tion of the infiaton potential from the parametrization of 
n s with n s — 1 = —a/N, and from the assumption be¬ 
tween the amplitude of the power spectrum — In A|j and 
n .j — 1 were also considered (Chiba 2015; Creminelli et al. 
2015; Gobbetti et al. 2015). On the other hand, the recon¬ 
struction of the infiaton potential from the power spectra 
by the method of functional reconstruction were usually 
applied (Hodges & Blumenthal 1990; Copeland et al. 1993; 
Liddle & Turner 1994; Lidsey et al. 1997; Peiris & Easther 
2006; Norena et al. 2012; Choudhury & Mazumdar 2014a; 
Ma & Wang 2014; Myrzakulov et al. 2015; Choudhury 
2016). 

In this paper, we discuss the reconstruction of the in¬ 
fiaton potential from the parametrizations of e(N), <fi(N) 
and n s (N). The parameters in the models are fitted to the 
Planck temperature and polarization data. The paper is or¬ 
ganized as follows. In the Section 2, we review the general 
relations between e(N), </(N), n s (N) and V(N) by apply¬ 
ing the slow-roll formula. The reconstruction of the infiaton 
potential from n s (N) is presented in Section 3. The recon¬ 
struction of the infiaton potential from e(N) is presented in 
Section 4. The reconstruction of the infiaton potential from 
</(N) is presented in Section 5. The conclusions are drawn 
in Section 6 . 


2 GENERAL RELATIONS 


For the single slow-roll inflation, to the first order of approx¬ 
imation, we have 


n a — 1 ~ 2r/ — 6e, r = 16e. 


( 1 ) 


Since 


dine 

dN 


2 V ~ 4e, 


SO 


n s 


-1 = 


— 2 e + 


dine 
dN ' 


( 2 ) 

(3) 


If we parametrize the slow-roll parameter e as a function 
of the number of e-folds before the end of inflation, then we 
can derive the parametrization n s (N) and r(N). Conversely, 
if we parametrize n 3 (N), we can solve equation (3) to get 
e(N), and then r(N). 

From the energy conservation of the scalar field, we have 


dlnp^, + 3(1 + Medina = 0. 


(4) 


Since p$ ~ V(</>), din a = —dN, and 


1 + 


£ 

P<t> 


<j> 2 2 

VW “ 3 6 


(5) 


1 d In V 

2 dN 


^(ln V),n > 0, 


( 6 ) 


and the inflton always rolls down the potential during infla¬ 
tion. Substituting the above result into equation (3), we get 
(Chiba 2015) 


ns “ 1 * "(In + ( ln ¥) , N = ( ln W) ■ (7) 

If we have any one of the functions e(N), n 3 (N) and V(N), 
we can derive the other functions by using equations. (3), 
( 6 ) and (7). 

To derive the form of the potential we need to 

find the functional relationship <f>(N) for the scalar field. 
Note that 


d(j> = ±y/27(N)dN, 


( 8 ) 


where the sign ± depend on the sign of the first derivative 
of the potential and the scalar field is normalized by the 
Planck mass M v i = (8-kG)^ 1 ^ 2 = 1, so 

r-N 

<t> - </>e = ± / s/2 7{N)dN. (9) 

Jo 

Once one of the functions e(N), n s (N), (/)(N) and V(N) is 
known, in principle we can derive n 3 , r and the potential 
V(4>) by using the relations (3), ( 6 ), (7) and ( 8 ). 

Before we present particular parametrizations, we 
briefly discuss the effect of second-order corrections. To the 
second-order of approximation, we have (Stewart & Lyth 
1993; Schwarz et al. 2001) 


n 3 — 1 « 2rj — 6 e — ( - 77 - + 24C ) L + ^77 


+ (16G - 2)ep + ( - - 2G ) C, 


16e 


1+ 4C-- e+ --2C r, 


^ * -16«7 + 24e 2 + 2£, 

where the Euler constant 7 ~ 0.577, C = ln 2+7 —2 
and 


( 10 ) 


( 11 ) 

( 12 ) 

-0.73 


(13) 


The observational data requires n B — 1 ~ 2 /N, so e and 77 are 
at most in the order of 1/N, their derivatives with respect 
to N have at most the order of 1/N 2 , so £ is at most in the 
order of 1/N 2 . Therefore, the running of the scalar spectral 
index and the second-order corrections will be in the order 
of 1/N 2 , we may neglect the second-order corrections. 


3 THE PARAMETRIZATION OF THE 
SPECTRAL TILT 


We approximate n 3 as 


n 3 — 1 fts — 


N + a ’ 


(14) 


where the constants p and a are both positive, and the con¬ 
stant a accounts for the contribution from the scalar field 
<j) e at the end of inflation. With this approximation, then n 3 
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is well behaved at the end of inflation when N = 0. Note 
that the functional form of the potential is not affected by 
a. The observational results favour p > 1, so we consider 
p > 1 only. 

For p > 1, the solution to Equation (3) is 


e(N) 


_ P- 1_ 

2(N + a) + C(N + ct)p’ 


(15) 


where C > 0 is an integration constant. This is a generaliza¬ 
tion of the Mukhanov parametrization e(N) = /3/(N + l) p 
(Mukhanov 2013). Since e(N = 0) tts 1 , so 2 a + Ca p « p— 1 . 
If C < 0, then |C| must be a very small number to ensure 
that 0 < e(N = 60) 1, and its contribution is negligible so 

that we can take it to be zero. Therefore, we consider C > 0 
only. Since 


C 



(16) 


so C > 0 requires that p > 2 a + 1 . 

Combining Equations (2) and (15), we get 


n(N) = -———--- 

/v ’ 2(N + a) + C(N + a)p 2 (N + a) 


(17) 


For N 1 and C ~ 1, e ~ l/(lV-|-a) p and p ~ — l/(7V-|-a), 
so the tensor to scalar ratio r ~ 1/(JV + oi) p is small and 
only p contributes to n s . 

Either solving Equation (7) with the parametrization 
(14), or solving Equation ( 6 ) with the solution (15), we get 
(Chiba 2015) 


V{N) 


p- i r i c 

A [ (N + a ) p_1 2 


(18) 


where A > 0 is an integration constant. 

Let us consider the special case C = 0 first. For C = 0, 
we get 


e(N) = 

P~ 1 

(19) 

2(N + a) ’ 

and 



p(N) = 

2p — 3 

( 20 ) 

2(N + a)' 


Therefore, both e and p contribute to the scalar spectral tilt 
and 


r = 


8 (P - 1) 
P 


(1 - n s ). 


( 21 ) 


Substituting Equation (19) into Equation (9), we get 
0 — 0 e = ±2\Jp — 1 (VN + a — y/a), (22) 



Planck TT.TE.EE+lowP 

- C=0,N-50 ‘ 

- C=0,N=60 

- p=2,N=50 

- P-2.N-60 

- 0(N),N-5O - 

0(N),K-6O - 


Figure 1. The marginalized 68 %, 95% and 99.8% CL contours for 
n s and ro .002 from Planck 2015 data (Planck Collaboration XX 
2015) and the observational constraints on the parametrizations 
(15) and = a In (3IV + 7 ). The cyan regions are for the 

parametrization (15) with C = 0, the brown regions are for the 
parametrization (15) with p = 2, and the magenta areas are for 
the parametrization = a In (/TV + 7 ). 


For simplicity, we apply the Planck 2015 68 % CL con¬ 
straints on r and n s to the r — n s relation ( 21 ) 
(Planck Collaboration XX 2015), and we find that no p sat¬ 
isfies the 68 % CL constraints as shown in Fig. 1. At about 
the 99.8% CL level, 1.00 < p < 1.93 and 0 < A0 < 12.39 
for N = 50, 1.20 < p < 2.21 and 6.65 < A0 < 15.42 for 
N = 60. From the above analysis, we see that the power-law 
potential is disfavoured by the observations at the 68 % CL, 
and the filed excursion for the inflaton is super-Planckian. 
For N = 50, p ~ 1 is marginally consistent with the ob¬ 
servational constraint at the 99.8% CL, it is possible that 
the field excursion of the inflaton is sub-Planckian and the 
tensor to scalar ratio r is close to zero. 

Next we consider the special case with p = 2. For this 
case, we get 


2y/2 


4>{N) = ± _ arcsinh\/C(Af + a)/2 + cj> 0 , 

vC 

2y/2 

4> e = db arcsinh^/Ca/2 + <j> 0 , 

VC 


(25) 

(26) 


where (j >0 is an arbitrary integration constant and C = (1 — 
2a)/a 2 > 0. Combining Equations (25) and (18), we get the 
corresponding T-model potential (Kallosh & Linde 2013), 


V{<j>) = tanh 2 


Vc_ 

2V2 


(0 - 4>o) 


= Vo tanh 2 [ 7 (<(> — <jf> 0 )], 


(27) 


or 

<P(N) = ±2V(P~1)(N + a) + 0o, (23) 


where (j >0 is an arbitrary integration constant, and ij> e = 
± 2fl(p — l)aA4>o- Since 2a ~ p— 1, so 1 and this model 
corresponds to large field inflation. Combining Equations 
(23) and (18), we get the power-law potential for chaotic 
inflation (Linde 1983) 


V{4>) = 


p- 1 


A(p - l) p - 


0—00 


2(p—1) 


= Vo(0- 0o) 2(p_1) . 


(24) 


where 7 = If C < {N + a) 1 or a —> 1 / 2 , then 

the potential reduces to the quadratic potential. From the 
discussion on the power-law potential, we know that the in¬ 
flaton is a large field and it is disfavoured by the observations 
at the 68 % CL level. If C > 1 or a 1, then the potential 
becomes 

V(0) « V&{l-2exp [- X /C72(0-0 O )]} 2 

« ■%, {1 - 4 exp [- y/cj2 (0 - 0o)] } . 

The potential includes the models with a-attractors 
(Kallosh et al. 2013) and the Starobinsky model 
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(Starobinsky 1980) when C = 4/3 or a = (y/21 — 3)/4. 
By fitting Equations (14) and (15) to the Planck 2015 
data, for N = 50, we get a < 0.4993 and 0 < A <p < 13.90 
at the 99.8% CL, for N = 60, we get a < 0.5009 and 
0 < A cp < 15.49 at the 99.8% CL. The results are shown 
in Fig. 1. Here we extend the integration constant C to 
the region of C < 0, and we verify the conclusion that 
|(7| is very small if C < 0 as discussed above. The above 
results also tell us that this model can be either small field 
inflation or large field inflation depending on the value of 
a. If a <C 1, then C > 1 and A (p is small. If a is close to 
1/2, then C is close to zero, and A (j> is large. 

For the general case with C > 0 and p / 2, we get 


4>(N) = h ± + q ) 1_P/2>< 

[1 p — 2 4 — 3 p 2{N + a) 1_p 

2 1 [ 2 ’ 2{p — 1) ’ 2 — 2p’ C 


The analytical form of the potential is not apparent, so we 
analyse the asymptotic form of the potential. For C C 1, 
the potential will be the same as the case with (7 = 0, and 
it is the power-law potential. For C > 1, Equation (15) can 
be approximated as 


e(A) 


P~ 1 

C{N + a) p ’ 


(30) 


and Equation (29) can be approximated as 


4 THE PARAMETRIZATION OF e(N) 


For the parametrization 
, . a 

1 + sexp(— /3N) ’ 

we get 


—2a + /3s exp (—/3N) 
1 + s exp(— /3N) 


( 34 ) 


(35) 


Note that for the parametrization in this section, we take 
s = ±1 for simplicity, a > 0, /3 > 0 and the parameter N = 
AL + N e so that we consider the contribution from <p e , i.e., at 
the end of inflation, 0 < N = N e -C A*. When s = — 1, we 
require exp(/3N e ) > 1. Substitute the parametrization (34) 
into Equation (6), we get 

V(N) = Ho[s + exp(/3 N)] 2a/fi . (36) 


Using Equation (8), we get 

</>(N) = ± 2v ^ In |cxp(/31V/2) + s + exp(/37V)j , 


and 

<t>e = (/)0 ± In ^exp(/3N e /2) + a/s + exp(/3AL)j . 

Combining Equations (36) and (37), we get 


%(/>) = Vo 


s + 




21 2 “/^ 


(37) 


(38) 


(39) 


<KN) = <Po± ^- \I 2{P , 1)C ? {N + a/ 2 -*^ 2 , (31) 

2 — p y p — 1 — 2a 


and 


= <t> 0 ± 


2a 2(p — 1) 


(32) 


2 — p | p - 1 — 2a 
Combining Equations (18) and (31), for (7 > 1, the potential 


V(<P) = 


2(p-l)a p 


A{p — 1 — 2a) 
2 — p 


it 


2 a p 


p — 1 — 2a 


1/(2 —p) 


1 —2(p-l)/(2—p) 


[2^=1 


{<t> ~ M 


(33) 


If p > 2, then for small held cp, the potential reduces to the 
hilltop potential V((p) = Vo[l — {<p/M) n \ (Boubekeur & Lyth 
2005) with n = 2(p — 1 )/{p — 2) (Garcia-Bellido & Roest 
2014; Creminelli et al. 2015). If 1 < p < 2, then for large 
held <p, the potential reduces to the form = Vo[l — 

(M/(p) n ] with n = 2(p— l)/(2 —p) (Garcia-Bellido & Roest 
2014; Creminelli et al. 2015). Fitting the model with general 
p and a to the observational data (Planck Collaboration XX 
2015), we hnd the constraints on p and a for N = 60 and 
the results are shown in Fig. 2. The results tell us that the 
model can accommodate both small and large held inhation. 
Because p > 2aTl, so the parametrization (14) requires that 
(l-n s )(N+a) > 2a+l. At the 99.8% CL, 0.02 <n s < 0.05, 
so a < 1.03 and 1.18 < p < 3.11 if we take N = 60. 


where U(4>) = exp(±/3(^> — </o)/(2\/2a)). For s = 1, the 
potential is 


V{<j>) = Vo 


{ /3((p - M 

V 2s/2a 


4ot/0 


(40) 


By htting Equations (34) and (35) with s = 1 to the Planck 
2015 data (Planck Collaboration XX 2015), we hnd that no 
a and j3 satishes the 99.8% CL constraints. 

For s = — 1, the potential is 


%(/») = Vo 


sinh 


{ P(<t> ~ <Po) \ 

V 2V2a J 


4 a/ 0 


(41) 


Fitting Equations (34) and (35) with s = — 1 to the Planck 
2015 data (Planck Collaboration XX 2015), we hnd that no 
a and /3 satishes the 68% CL constraints, so the model is dis¬ 
favoured at the 68% CL. The 95% and 99.8% CL constraints 
on a and /3 for N = 60 are shown in Fig. 3. By using the 
99.8% CL constraints, we hnd that 8.37 < Arp < 16.48, so 
the held excursion of the inhaton in this model is super- 
Planckian. 

For the parametrization 


e(N) 


a exp(— /3N) 

1 + sexp(—f3N) ’ 


(42) 


we get 

_ /3 T 2aexp(— /3N) 

ns 1 + s exp (—/3N) 


(43) 


Substitute the parametrization (42) into Equation (6), we 
get 

V(N) = V 0 [l + sexp(-/3N)]- 2a/(f,s) . (44) 
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n s a 

Figure 2. The marginalized 68%, 95% and 99.8% CL contours for n s and ro .002 from Planck 2015 data (Planck Collaboration XX 2015) 
and the observational constraints on the parametrization (15). The left-hand panel shows the n s — r contours and the right-hand panel 
shows the constraints on p and a for N = 60. The green, yellow and red regions correspond to 68%, 95% and 99.8% CLs, respectively. 
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Figure 3. The marginalized 68%, 95% and 99.8% CL contours for n s and ro .002 from Planck 2015 data (Planck Collaboration XX 
2015) and the observational constraints on the parametrizations (34) with s = — 1. The left-hand panel shows the n s — r contours and 
the right-hand panel shows the constraints on a and (3 for IV = 60. The yellow and red regions correspond to 95% and 99.8% CLs, 
respectively. 



From Equation (8), we get 


,, N , 2y/2 a 

<p{N) = 0o =F--=-arcoth 

PVs 

For s = 1, we get 

2fl2a 


s + exp (/37V) 


= <Po T ■ 


ft 


arcoth ^ \/l + exp(/3N e ) j 


and the potential 


V(4>) = Vb 


sech 


2V2a 


4a/(3 


(45) 


(46) 


(47) 


If /3 = 4a = 2/M 2 , we recover the potential for the s-dual in¬ 
flation V = Vosech (Anchordoqui et al. 2014). Fitting 
Equations (42) and (43) with s = 1 to the Planck 2015 data 
(Planck Collaboration XX 2015), we obtain the constraints 
on the parameters a and /3 for IV = 60 and the results are 
shown in Fig. 4. From Fig. 4, we see that the s-dual infla¬ 
tion is consistent with the observational data. By using the 
99.8% CL constraints, we find that 0 < A</> < 26.57, so the 


model includes both the large field and small field inflation. 
If a is close to zero, then A <j> is small. 

For s = — 1, we get the potential, 


V{<j>) = Vo 


V 2^2 a 


4a/ (3 


(48) 


or 


v(4>) = Vo 


Vo 

4 


( P{<t>- <M 

V 2\[2a 


4 a//3 


i + cos 

V V?R 


2a//3 


(49) 


If we take /3 = 2a and / = l/\/2a, then we recover the 
potential for natural inflation V(0) = A 4 [l + cos(</>//)] 
(Freese et al. 1990). Fitting Equations (42) and (43) with 
s = —1 to the Planck 2015 data (Planck Collaboration XX 
2015), we obtain the constraints on the parameters a and 
/? for N = 60 and the results are shown in Fig. 5. From 
Fig. 5, we see that natural inflation is disfavoured at the 
68% CL. By using the 99.8% CL constraints, we find that 
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Figure 4. The marginalized 68 %, 95% and 99.8% CL contours for n s and ro .002 from Planck 2015 data (Planck Collaboration XX 
2015) and the observational constraints on the parametrizations (42) with s = 1. The left-hand panel shows the n s — r contours and the 
right-hand panel shows the constraints on a and (3 for N = 60. The green, yellow and red regions correspond to 68 %, 95% and 99.8% 
CLs, respectively. The solid line (3 = 4a in the right-hand panel corresponds to the s-dual inflation. 


0 < A(j) < 54.41, so the model includes both the large 
field and small field inflation, and the small field inflation 
is achieved when a is close to zero. 

For the parametrization 


e(N) = aex P(~P N ) 

1 ’ [l + sexp(—/3A0] 2 ’ 

we get 

—ft — 2aexp(—/37V) + /3s 2 exp(—2/3iV) 
[1 + sexp(—ft IV )] 2 


(50) 


(51) 


_2rv 

V (N) = Vo exp 
and 

<KN) = <fto ± r 

pV s 

For s = — 1, we get the potential 


_ft[s + exp(ftiV)] J ’ 

2\f2a arc t an [exp(/?A r /2)/v / s] • 


V ((ft) = Vo exp 


-ggsinh 


ft 


V 2y/2a 


(52) 


(53) 


(54) 


If ft = 2a = 4/p 2 , then the potential reduces to the hilltop 
potential V(<ft) fa Vo [1 — (4‘/p) P ] with p = 2 for small <j>. If 
a = /3 = 8/p 2 , then the potential reduces to the double well 
potential V(cp) = Vo[l — (0//z) 2 ] 2 for small cj> (Olive 1990). 
Fitting Equations (50) and (51) with s = —1 to the Planck 
2015 data (Planck Collaboration XX 2015), we obtain the 
constraints on the parameters a and j3 for IV = 60 and the 
results are shown in Fig. 6 . From Fig. 6 , we see that the 
double well potential is excluded by the observational data 
and the hilltop potential with p = 2 is disfavoured at the 
68 % CL. By using the 99.8% CL constraints, we find that 
0 < A (f> < 14.31, so the model includes both the large field 
and small field inflation. If a is close to zero, then A<j> is 
small. 

For s = 1, we get the potential 


V(</>) 


Vo exp 



^ 1 + cos 


ft Oft - 0o) V 
)_ ' 


(55) 


Fitting Equations (50) and (51) with s = 1 to the Planck 


2015 data (Planck Collaboration XX 2015), we obtain the 
constraints on the parameters a and ft for N = 60 and 
the results are shown in Fig. 7. By using the 99.8% CL 
constraints, we find that 0 < A (ft < 31.23, so the model 
includes both the large field and small field inflation, and 
the small field inflation is achieved if a is close to zero. 


5 THE PARAMETRIZATION OF <ft(7V) 

Combining Equations ( 6 ) and ( 8 ), we get 

(In V), N = (</>, N ) 2 . (56) 

Once the functional form <ft(lV) is known, we can derive the 
potential form E(<ft). Let us first consider the power-law 
parametrization 

<t>(N) = a(X + 7 ) /3 . (57) 

For ft = 1/2, from Equation (56), we get the power-law 
potential, 

ct 2 

V(<ft) = Vo(N + 7 )T = Eo 2 , (58) 


where Vo is and integration constant. From Equation ( 8 ), we 
get 


8(N + 7 )’ 

so 

1 + ct 2 /4 2 a 2 

ns ~ ~ A + 7 5 r N + 7 ’ 


(59) 


(60) 


and a 2 m 87 . The results are the same as those discussed in 
Section 3 with C = 0, p = 1 + a 2 /4 and a = 7 . 

For ft 5 ^ 1/2, we derive the potential, 


V (ift) = Vo exp 
= Vo exp 


' U~ 2 ft 2 
2 ft — 1 


(N + 7 ) 2/3 " 1 


a 2 p 2 

2ft -1 W 


(61) 
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a 


Figure 5. The marginalized 68%, 95% and 99.8% CL contours for n s and ro .002 from Planck 2015 data (Planck Collaboration XX 2015) 
and the observational constraints on the parametrizations (42) with s = — 1. The left-hand panel shows the n s — r contours and the 
right-hand panel shows the constraints on a and (3 for N = 60. The green, yellow and red regions correspond to 68%, 95% and 99.8% 
CLs, respectively. The solid line (3 = 2a in the right-hand panel corresponds to the natural inflation. 



n s a 

Figure 6. The marginalized 68%, 95% and 99.8% CL contours for n s and ro .002 from Planck 2015 data (Planck Collaboration XX 2015) 
and the observational constraints on the parametrizations (50) with s = — 1. The left-hand panel shows the n s — r contours and the 
right-hand panel shows the constraints on a and (3 for N = 60. The green, yellow and red regions correspond to 68%, 95% and 99.8% 
CLs, respectively. The black solid line /3 = 2a in the right-hand panel corresponds to the hilltop potential with p = 2, and the magenta 
solid line (3 = a corresponds to the double well potential. 


the scalar spectral tilt, 


n s l- Ar + 7 oft{N + 7 ) p , 

(62) 

the tensor to scalar ratio, 


r = 8a 2 (3 2 (N + 7 ) 2 ^ -2 , 

(63) 


and the parameters cr, (3 and 7 satisfy the relation a 2 /3 2 ~ 
2 / y 2 ~ 2/3 . Note that Equations (62) and (63) include the spe¬ 
cial case (3 = 1/2 which corresponds to the power-law po¬ 
tential. Fitting Equations (62) and (63) for the parametriza- 
tion (57) to the Planck 2015 data (Planck Collaboration XX 
2015), we obtain the constraints on the parameters a and (3 
for N = 60 and the results are shown in Fig. 8 . By using the 
99.8% CL constraints, we find that 0 < A</> < 754.23. If (3 
is near zero, then it is large field inflation. If 7 is near zero, 
then it is small field inflation. In the left-hand panel of Fig. 
8 , the lower bounds of the n s — r contours are set by (3 = 0. 

Next we consider the parametrization (f>(N) = 


a ln(/37V + 7 ). Combining Equations (1), (3), ( 8 ) and (56), 
we get 

V(fl = V 0 exp(-^A_) 

= Vo exp [-a 2 /3exp(—0/cr)] , 

_2 q2 

/ ht\ cr p 

e( ) = 2(/3N + 7) 2 ’ 


s ~ N + 7 //3 (7V + 7//3) 2 ’ 

8 cr 2 

r " (N + 'y/ft) 2 ■ 

From the end of inflation condition e(N = 0) « 1, we get 
cr 2 /? 2 = 2y 2 , so r = 16(7 //3) 2 /(N + 'y//3) 2 , and n s and r 
depend on the parameter 7 //? only. If 7 /ft is small, then 
n 3 — 1 ~ —2/N and r ~ 1 /N 2 , and the model will behave 
like the model (14) with p = 2 and small a = 7 //3, so the 


(64) 

(65) 

( 66 ) 
(67) 
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Figure 7. The marginalized 68 %, 95% and 99.8% CL contours for n s and ro .002 from Planck 2015 data (Planck Collaboration XX 
2015) and the observational constraints on the parametrizations (50) with s = 1. The left-hand panel shows the n s — r contours and the 
right-hand panel shows the constraints on a and ft for N = 60. The green, yellow and red regions correspond to 68 %, 95% and 99.8% 
CLs, respectively. 




i 




n s y 

Figure 8. The marginalized 68 %, 95% and 99.8% CL contours for n s and ro .002 from Planck 2015 data (Planck Collaboration XX 2015) 
and the observational constraints on the parametrizations (57). The left-hand panel shows the n s — r contours and the right-hand panel 
shows the constraints on (3 and 7 for N = 60. The green, yellow and red regions correspond to 68 %, 95% and 99.8% CLs, respectively. 


two models will cover some common regions in the n s — r 
graph. The fitting results are shown in Fig. 1. In particular, 
for (j)/cr 1 , the potential can be approximated as 

V(<j>) « Vo [l - 2 y 2 exp(—0/cr)l 

, ( 68 ) 

« Vo [l - 7 exp (~(j)/a)\ . 

Therefore, the R 2 inflation is also included in this model. 

In the last we consider the exponential parametrization 
<j>(N) = aexp(/3N + 7 ) with /3 < 0. Following the same 
procedure as above, we obtain 

V{4>) = Vo exp(/3<)> 2 /2), (69) 

e = a ^ exp(2/3iV + 27 ) = r/16, (70) 

n s — 1 = 2/3 — /3 2 a 2 exp(2/3N + 27 ). (71) 

The model parameters satisfy the relation /3 2 cr 2 exp( 27 ) = 2, 
so e = exp(/31V) and n s — 1 = 2/3 — 2exp(2/3N). Since e < 1, 
the parameter /3 should be negative. From the constraint on 
r, we can get the upper limit on /3, with this upper limit, we 


find that the model is not consistent with the observational 
data. 


6 CONCLUSIONS 

For the double well potential V(< p) = Vo[l — {cj>/p) 2 ] 2 , 
the potential V(rf>) = Vo exp(/3(j> 2 ) and the potential 
V{4>) — Vo[cosh(/30/2v / 2Q)] 4 “^ /3 , the predicted n s and r 
are not consistent with the Planck 2015 data. The power- 
law potential V{(j>) = Vb0 2p-2 , the potential V{cj)) = 
Vb[sinh(/30/2-\/2a)] 4 “^ /3 , the natural inflation and the hill¬ 
top potential with n = 2 are disfavoured by the observa¬ 
tional data at the 68% CL. At the 99.8% CL, we find that 
1.2 < p < 2.21 for the power-law potential if we take the 
number of e-folds before the end of inflation N = 60. For 
the power-law potential V(4>) = Vo((p — the T- 

model potential V(4>) = Vo tanh 2 ( 7 ()>) which includes the 
a-attractors and the Starobinsky model, the hilltop poten¬ 
tial V(cf>) = Vo[l — (</>/M) n ] with n = 2(p — 1 )/{p — 2) 
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and p > 2, the potential V(0) = Vo[l — (M/0) 71 ] with 
n = 2 (p — 1 )/ (2 — p) and 1 < p < 2 , and the potential 
V( 0 ) = Vo exp[—<r 2 /3exp(— 0 /cr)], their spectral tilts have 
the universal behavior n s = —p/(N + a). 

For the parametrization e(N) = aexp(—/3N)/[l + 
exp(—/31V)], we get n s — 1 = —\/3 + 2a exp(— j3N)]/[l + 
exp(— /3N)\ and the corresponding potential V(<j >) = 

Vo[sech(/30/2\/2a)] 4c “ ,/ ' S which includes the s-dual infla¬ 
tion. For the parametrization e(A r ) = aexp(—/3N)/[l — 
exp(—/3N)], we get n s — 1 = — [/3 + 2a exp(— (3N)\/[1 — 
exp(— /3N)\ and the corresponding potential V(<j >) = 

Vb[sin(/30/2%/2a)] 4 “^ which includes the natural infla¬ 
tion. For the parametrization e(N) = aexp(—/3N)/[l — 
exp(—(3N)] 2 , we get the corresponding potential V(0) = 
Vo exp[—2a sinh 2 (/30/2\/2a)//3] which includes the hilltop 
potential with n = 2 and the double well potential. The 
tensor to scalar ratio r for these models can easily be small 
due to the factor exp(— /3N) in e(N). For the parametriza¬ 
tion 0(A) = a\n(/3N + 7 ), the corresponding potential is 
V(0) = Voexp[—CT 2 /3exp(—0/cr)], both n s and r depend 
on 7 //? only and the model has the universal behavior 
n s — 1 = —2/{N + 7 /ft) if 7 //? is small. All these models can 
achieve both small and large field inflation. 

Based on the slow-roll relations (3), ( 6 ), (7) and ( 8 ), 
by parameterizing one of the parameters n s (N), e(N) and 
0(A), and fitting the parameters in the models to the ob¬ 
servational data, we not only obtain the constraints on the 
parameters, but also easily reconstruct the classes of the 
inflationary models which include the chaotic inflation, T- 
model, hilltop inflation, s-dual inflation, natural inflation 
and R 2 inflation, and the reconstructed inflationary models 
are consistent with the observations. Since the observational 
data only probes a rather small intervals of scales, the re¬ 
constructed potentials approximate the inflationary poten¬ 
tial only in the slow-roll regime for the observational scales 
10 ~ 3 Mpc<fc - 1 <10 4 Mpc. Outside the slow-roll regime, 
the inflationary potential can be rather different, but it does 
not mean that the reconstructed potential is not applicable 
in that regime. Once the potential is obtained, we can either 
apply the slow-roll formulae or work out the exact solutions. 
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